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PREFACE 

The  work  presented  In  this  report  vas  done  under  NASA  Contract 
NASr-21(02),  monitored  by  the  Director  of  Communication  Studies,  in  the 
Office  of  Applications.  The  general  study  area  of  the  contract  comprises 
technical  studies  relating  to  communication  satellites.  As  part  of  this 
work,  RAND  has  paid  particular  attention  to  guidance  and  orbit  stabiliza¬ 
tion  problems  of  24-hour  synchronous  satellites.  The  present  study  shows 
that  the  problem  of  maintaining  such  satellites  In  cm  exactly  stationary 
position  will  be  more  complex  than  had  been  previously  expected.  The 
results  should  be  of  concern  to  all  agencies  and  contractors  involved  in 
the  development  of  a  synchronous  satellite. 

A  summary  of  the  findings  of  this  study  vas  presented  to  NASA  on 
October  20,  1961. 
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SUMMARY 

This  paper  presents  the  results  of  an  Investigation  of  the  behavior 
of  a  synchronous  (24 -hour)  satellite  as  effected  by  the  triaxiality  of  the 
earth.  This  includes  not  only  the  effect  of  the  equatorial  bulge  but  also 
the  ellipticity  of  the  earth's  equatorial  section.  The  results  indicate 
that  there  are  only  two  positions  on  the  equator  (123°J'  West  Longitude 
and  56°J1'  East  Longitude)  at  which  a  truly  synchronous  satellite  can  exist 
in  a  stable  condition.  In  order  to  establish  a  synchronous  satellite  at 
any  other  longitude,  station- keeping  propulsion  of  the  order  of  51  ft/ sec  for 
each  year  of  operation  would  be  required.  If  no  station-keeping  propulsion 
is  provided,  the  satellite  will  execute  long  period  (greater  than  1.3  years) 
oscillations  about  the  closest  of  the  two  stable  positions  mentioned  above. 
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I.  DURODUCTICH 

In  recent  years  the  Idea  of  establishing  an  artificial  satellite  in  a 
synchronous  equatorial  orbit  about  the  earth  has  become  increasingly  at* 
tractive.  Since  such  a  satellite  by  definition  would  remain  above  the  same 
point  on  the  earth's  equator ,  it  could  be  used  as  a  communication  relay  station 
between  any  two  points  on  the  earth  which  are  within  lta  field  of  view. 

If  it  is  assumed  that  the  mass  distribution  of  the  earth  is  spherically 
synmetric ,  then  the  resulting  gravitational  potential  varies  inversely  as  the 
radial  distance  from  the  earth's  center  and  is  independent  of  either  latitude 
or  longitude.  Under  this  assumption,  it  is  a  simple  matter  to  demonstrate 
that  a  satellite  on  a  circular  orbit  at  an  altitude  of  22,236  st  mi  would 
have  an  orbital  period  of  one  sidereal  day.  Thus,  such  a  satellite  in  the 
earth's  equatorial  plane  would  appear  to  be  fixed  relative  to  the  earth. 

It  has  long  been  recognized  that  in  order  to  take  into  account  the  actual 
non-spherlcal  nature  of  the  earth's  mass  distribution,  it  is  necessary  to 
Include  additional  terms  in  its  potential  function  in  the  form  of  higher 
order  spherical  harmonics.  These  additional  terms  depend  not  only  on  the 
radial  distance  but  also  on  latitude  and/or  longitude.  As  a  result  of  obser¬ 
vations  of  some  of  the  artificial  satellites  launched  thus  far , ^  it  has 
been  possible  to  refine  the  estimates  of  the  magnitude  of  these  additional 
spherical  harmonic  terms.  Ike  two  additional  terms  considered  in  the  above 
reference  are  (l)  the  zonal  solid  spherical  harmonic  corresponding  to  the 
usual  oblateness  of  the  eaurth  such  that  the  equatorial  radius  exceeds  the 
polar  radius,  and  (2)  the  sectorial  solid  spherical  harmonic  which  results 
from  the  ellipticity  of  the  eairth's  equatorial  section.  The  first  of  these 
terms  is  latitude  dependent  while  the  second  is  a  function  of  longitude  as 


well 
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The  question  now  arises  as  to  vhat  effect  these  modifications  in  the 
earth's  potential  might  have  on  the  synchronous  nature  of  a  24-hour  satellite. 
The  present  paper  addresses  this  problem  by  setting  up  the  general  equations 
of  motion  of  a  satellite  under  the  Influence  of  the  expanded  potential 
function,  These  equations  are  then  linearized  in  terms  of  small  deviations 
from  a  circular  synchronous  orbit.  On  the  assumption  that  the  orbital 
injection  procedure  is  perfect  and  that  the  initial  deviations  from  the 
synchronous  condition  are  all  zero,  the  above  equations  can  be  solved  for 
the  forced  solution  resulting  from  the  added  terms  in  the  potential  function 
which  act  as  driving  functions  in  the  linearized  equations. 

On  the  basis  of  these  results  it  is  possible  to  determine  the  amount 
of  propulsion  required  to  maintain  a  synchronous  satellite  at  a  given  longitude 
on  the  equator. 

While  the  solution  of  the  linearized  equations  is  only  valid  for 
variations  in  orbital  radius  and  for  variations  in  epoch  angle  up  to  10° 
relative  to  the  synchronous  position,  an  approximate  equation  is  also  develop¬ 
ed  which  allows  large  variations  in  epoch  angle  but  is  6till  restricted  to 
small  radial  displacements.  By  means  of  thiB  equation  the  behavior  of  the 
satellite  can  be  described  in  the  event  no  propulsion  is  provided  to  main¬ 
tain  position. 
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II.  METHOD  OF  A1AIY8I3 


STATEMENT  OF  THE  PROBLEM 

The  problem  which  is  solved  in  the  succeeding  sections  of  this  report 
is  the  determination  of  the  motion  of  a  presumably  synchronous  satellite 
relative  to  the  earth  when  the  effects  of  earth's  oblateness  and  the  elllp- 
tlcity  of  its  equatorial  section  are  Included  in  the  gravitational  potential 
function.  For  the  purposes  of  this  analysis ,  the  gravitational  attraction 
of  the  sun  and  moon,  the  sun's  radiation  pressure,  and  residual  drag  effects 
are  neglected. 

EQUATIONS  OF  MOTION 
Coordinate  System 

The  reference  system  adopted  is  shown  in  Fig.  1  where  the  XYZ  system  is 
an  inertial  reference  with  origin  at  the  earth's  center  and  OZ  along  the  polar 
axis.  In  this  system  the  position  of  the  satellite  P  is  specified  by  the 
spherical  coordinates  0,  0,  and  r.  Also,  the  line  OA  represents  the  instan¬ 
taneous  position  of  the  minor  axis  of  the  earth's  equatorial  section  so  that 
the  angle  y  represents  the  difference  in  longitude  between  the  satellite  and 
this  minor  axis  reference  line.  For  a  truly  synchronous  satellite  y  should 
remain  constant.  In  addition,  three  unit  vectors,  r^,  n^  and  n^  are  specified 
at  P  corresponding  to  displacements  due  to  small  changes  in  r,  0,  and  0,  re¬ 
spectively. 

Acceleration  Components 

The  components  of  acceleration  in  the  direction  of  the  three  unit  vectors 
specified  above  can  be  expressed  in  terms  of  r,  0,  and  0  as  follows 


•  •  •  o  2  *  2 

a  -  r  -  r®  cos  0  -  r0 
r 


a) 


Fig.  I  —  Coordinate  system 
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Force  Components 


The  formulation  of  the  earth's  gravitational  potential  presented  in 
Reference  1  is  used  to  determine  the  corresponding  force  components. 


goRo  r 

Uo  o  ■* 
r 


,  R0 

J2j 

r 


sin  0-1, 


|2)  **Q  O 

♦  3Jg  c08  0  008  27 


where 


R  ■  earth's  radius 
o 

g  »  gravitational  acceleration  at  the  earth's  surface 
J2  «  non-dimensional  multiplier 
J  (2) 

2  a  non-dimensional  multiplier 

In  Eq.  (U)  the  first  term  results  from  the  spherical  part  of  the  earth's 

mass  distribution,  the  second  or  term  from  the  usual  oblateness  about 

(2) 

the  polar  axis  and  the  third  or  Jg'  '  term  from  the  ellipticity  of  the 

earth's  equatorial  section.  It  is  seen  that  the  J2  term  depends  only  on 

(2) 

the  latitude,  0,  while  the  J2  term  depends  not  only  on  the  latitude  but 
also  on  the  relative  longitude,  7.  See  Fig.  2. 


The  numerical  values  of  J2  and  J2 


have  been  determined  in  Ref.  1  on 


the  basis  of  observations  of  Vanguards  II  and  III  as  follows 

J2  -  ♦  1,08219  x  10’3 


-  5.35  x  10 


t 
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A  value  of  J^2)  of  tiliB  magnitude  is  equivalent  to  a  difference  of  671  ft 
between  the  semi-major  and  semi-minor  axes  of  the  equatorial  section. 

The  desired  force  components  can  be  obtained  from  Eq.  (4)  as  follows 


,  *2  soRo 
- - 4— 


2r 


(3  »in20 


1) 


a  1  au 

b  *  ?3e 


9J2(2)  goR04  2 
— = - jp-£_  cos  0  cos  27 


oj  ^  g  R  4 

- - -  cos2*  sin  27 


(5) 

(6) 


fe  m  1  <Xj 

m  "  r  3* 


*2  goRo 


sin  0  cos  0 


where 

m 

7 


“T 

r 


cos  0  sin  0  cos  27 


mass  of  the  satellite 
6  ’  ®E 


(7) 


Complete  Equations 

The  general  form  of  the  equations  of  motion  is  obtained  by  equating 
the  corresponding  acceleration  and  force  expressions  from  E^s.  (l),  (2),  (3), 
(r0>  (u)  and  (7)  to  give  the  following 


'*  *2  2  *2 
r  -  r©  cos  0  -  r* 


o  o 
2 


o  o 


2r 


(3  sin  0-1) 


6 


9j2(2>  g0R04  2 

L' '  COB  0  COB  27 


(8) 


1  d  r_2l  _ 2^1  ^2  6 

ri3T5  It  lr  6  CM  *\  * - T» 


6j  ti  ‘  , 

2  °  °-  COB2#  Bln  27 


(9) 


Id  f  2*  1  *2  3J« 

r  dtLr  +  ^  cob  #  Bin  #  •  -  ■  Bln  0  cob  # 


Wa) 


sin  #  cob  0  cos  27 


(10) 


Since  the  present  analysis  is  concerned  with  an  equatorial  satellite,  the 
above  equations  can  be  simplified  by  assuming  that  0  is  a  snail  angle,  in 
which  case  Eqs.  (£),  (9)  and  (10)  reduce  to 


ae  *p 

r  -  re  -  - 


goRo2  3J2  8 A4  9J2(2)  «oRo4 

~5 - - -5 -  eo.  2, 

r  2r  r 


(11) 


1 
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1 
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d 

fr2al 

*2  *0*0 

dt 

L 

“  "  k 

r 

d 

r  2*  i 

*2 

dt 

[r  0j 

+  r0  0  ■  0 

sin  27 


(12) 

(13) 


( 2 ) 

where  terms  of  the  order  of  Jg#  and  JgV  '0  have  been  neglected  in  Eq.  (13). 

(2) 

Since  the  Jg  and  Jg  terms  do  not  have  any  significant  effect  on 
the  0  equation,  the  rest  of  the  analysis  will  be  concerned  only  with  notion 
in  the  orbital  plane,  as  described  by  Eqs.  (11)  and  (12). 


LINEARIZED  EQUATIONS 

Since  this  paper  is  primarily  concerned  with  deviations  of  the  satel* 
lite  from  a  nominally  synchronous  orbit,  it  is  convenient  to  linearize 


Eqs.  (11)  and  (12)  in  the  vicinity  of  the  conditions  for  such  an  orbit 
This  is  accoapliahed  by  defining  the  following  set  of  variables 
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r  -  rc  *  dr  (14) 

B  *  ®E  +  7o  *  Ly  (*5) 

7  •  7q  *  L7  (16) 


where 

rc  ■  orbital  radius  of  a  synchronous  orbit 

7q  •  longitude  difference  between  the  minor  axis  of  earth's 
equatorial  section  and  the  desired  synchronous  position 
of  the  satellite. 


Substitution  of  Eqs.  (14),  (15)  and  (16)  into  Ec.  (11)  gives 

2 


Ar  *  rceE  *  ®e  4r  -  ®e  •  *  -fr  (1  •  fr) 

c  c 


gR 
o  o 


2dr, 


3J2  6oRo  *+drx 

- jj -  U  -  — ) 

2  r  rc 


9J2(2)  g«R,  4 


^  (1  -  f 

r  4  rc  L 


cos  27  -  sin 

o 


2  ] 


(17) 


By  equating  the  steady  state  terms  to  zero  as  follows 


Ve  -  G°R°2 


3J2  «qRq 

2  rk 


=  c 


(li.) 


a  relationship  is  obtained  which  defines  the  value  of  the  synchronous  orbital 

radius,  r  .  Thus  it  is  seen  that  earth  oblateness  modifies  the  value  of 
'  c 
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r  from  that  obtained  for  a  ipberlcal  earth, 
c 

By  combining  Eqs.  (17)  and  (18)  and  neglecting  higher  order  terms 
the  following  equation  results 

T-o  *  2 

dr 


r-o  2  g  R  -| 

-  ♦  -^rj  *r  - 2  rc  ®i  47 


.  <2)  u 

1W2  8oR0  ill  ar0 


dr 


(2) 

9J2  *o"o 


COS  27 


(19) 


C  C 

Equation  (19)  can  be  further  simplified  by  Introducing  the  following  lumped 


constants 


(20) 


k2  <?  2 

rc 


(a) 


and  by  making  use  of  the  approximate  form  of  Eq.  (18)  in  which  the  term 
Is  neglected  to  give 

gJL2 

(22) 


•9 

O  0 


Thus  Eq.  (19)  becomes 


r  •  3»e  r  •  »e  47  *  l6ll2  ®i  •1"  % 47  - 

c  c 

+  9k2  9*  COB  27q 

Finally,  if  a  non-dimensional  time  t  Is  defined  by  the  relation 


(23) 


t  -  9£t 


(24) 
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Kq.  (23)  can  be  rewritten  as 

^  (f)  -  3(^)  -  2  *  (I8k2  .in 

dT  c  c 

In  a  completely  analogous  manner  Eq.  (12) 
non-dimensional  form  as  follows 

2  h  <r>  *  (3°k2  •*“  2V  r 

c  c 

♦  ^5  -  (12k-  cos  27  )  A7 


27q)  A7  ■  ♦  9fcg  cos  27q  (25) 

can  also  be  linearised  and  put  into 


-  *  6k0  sin  27 
2  o 


(26) 


The  preceding  analysis  has  reduced  the  two  non-linear  differential  equations 
of  motion  (11)  and  (12)  to  two  coupled  linear  differential  equations  (25) 
and  (2o)  with  constant  coefficients  relating  the  variables  Ar/rc  and  A 7. 

By  combining  Eos.  (25)  and  (26)  separate  equations  for  Ar/rc  and  A7 
are  obtained  in  operator  form  as  follows 


s4  ♦  s2  ♦  (2-*k2  sin  270)s  ♦  (j6k2  cos  27q)  ^  -  -  106k|  (27) 

J  c 

S4  +  s2  +  (24k2  sin  27q)  s  «■  (36kg  COS  27c)j  A7  -  -  16kg  sin  27q  (26) 


where  s 


d_ 

dt 


It  is  interesting  to  note  at  this  point  that  the  above  equations  depend 

only  upon  k2,  which  is  related  to  the  ellipticity  of  the  equatorial  section. 
The  only  significant  effect  of  the  other  oblateness  term  is  to  modify  the 
value  of  rc  through  Eq.  (16). 

In  order  to  solve  Eqs.  (27)  and  (26)  the  roots  of  the  characteristic 
operator  equation  are  obtained  in  the  form 
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•x  ■  ♦  12kg  tin  2 70  ♦  J  (29) 

•g  .  ♦  12kg  .in  27q  -  J  (30) 

#3  -  -  12kg  .In  2 70  ♦  6y."kg  cos  27Q  (31) 

%  -  -  12k2  iin  “  6<v/‘k2  co*  27 o  (32) 

The  fora  of  the  solution  depends  on  the  sign  of  cos  2y  since  this 

o 

determines  whether  and  are  real  or  conjugate  caaplex.  If  cos  2?o<0 
then  the  solution  is  of  the  form 

41  -  »2  [l  *  7  .1»2  2r0]  V2  .  ♦  (ia3  ,ln  V’  CO.  (,-«,) 


-  2  tan  2>q  >/-kg  cos  27^  e 


-  (I2k_  sin  27. )t 


*2 


cos  27 


*2  "'o'  sinh  (6y-kg  cos  27q)t 

(33) 


A7  •  -  6kg  [l  +  7  sin2  270]  ^  e  *  ^Uk2  8ia  ^o^  sin  (t-c^ 

*2  ”***  "'o)T  cosh  (6y-kg  cos  27q)t 


1  .  -  (l2k„  sin  27 Jt 


♦  ^  tan  27o  e 


-  ^  tan  27c 


(3M 


where 


tan  «i  -  3  tan  270 


(35) 


In  obtaining  the  above  solution,  it  has  been  assumed  that  the  initial  orbital 

•  • 

injection  errors  have  been  reduced  to  zero  so  that  (hr  *  hr  *  by  *  by  *  o) 

o  o  o  9  o 

and,  as  a  result,  Eqs.  (33)  and  (34)  represent  the  driven  solution  resulting 

(2) 

from  the  J,/  term  in  the  earth's  potential  function. 
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An  examination  of  the  above  solution  shove  that  a  number  of  simplifica¬ 
tions  can  be  made  in  view  of  the  magnitude  of  the  constant  k2  (of  the  order 
of  10“7).  "Hie  exponential  terms  have  time  constants  in  excess  of  300  yr 
and  can  be  replaced  by  unity.  Likewise,  terms  involving  kg  as  a  multiplier 
can  be  neglected  compared  to  those  involving  kg^2.  Thus,  a  simplified 
version  of  the  solution  can  be  written  as 

■  -  2  tan  2>q  k2  cos  2yQ  sinh  (6  ^/-k2  cos’"??  )t  (36) 

c 

A7  =  -  |  tan  2yQ  -  cosh  (o  •/-  k?  cos  270)tj  (37) 

This  form  of  the  solution  applies  for  *+5°  <|ro|<9°0. 

If  the  value  of  cos  27  q>  0,  then  the  roots  s^  and  3^  become  complex 
and  it  can  be  shown  that  the  form  of  the  solution  analogous  to  Eqs.  (36)  and 
(37)  is  given  by 

~  »  2  tan  27£jv/k2  cos  2?q  sin  (q^/kg  cos  27q)t  (36) 

8 

A7  =  -  \  tan  27q  -  cos  (cy/kg  cos  27q)tJ  (39) 

This  form  of  the  solution  applies  for  |7Q|  <  *+5°. 

Due  to  the  linearization  procedure  used  in  obtaining  the  above  solutions, 
they  are  only  valid  as  long  as 

dr«r 

c 

A  7  f  sin  A7 

An  examination  of  Eqs.  (30)  and  (36)  indicates  that  the  first  condition 

will  be  satisfied  for  values  of  t  equivalent  to  a  three  year  duration,  at 

which  time  dr  would  be  about  one  per  cent  of  r  . 

c 


However,  the  buildup  of  by  la  sore  rapid  and  nay  reach  a  value  of  10° 
In  about  two  months .  Thus  it  la  seen  that  the  behavior  of  by  la  the  de¬ 


termining  factor  with  regard  to  the  time  interval  over  which  the  aolutloca 
remain  valid. 

LARCH  AMPLE  EQUATIONS 


In  view  of  the  fact  that  by  exceeds  the  linear  range  In  the  above 
solution  after  a  reasonable  length  of  time,  it  was  felt  that  a  solution 
which  was  valid  for  large  values  of  by  would  be  desirable.  This  can  be 
obtained  from  Eqs.  (11)  and  (12)  by  substituting  +  y  for  9  to  give 


•  •  •  •  o  8 

r  -  r  (S  .  ,)2  -  -  -y- 


»2  «oBo 


cos  27  (  40) 


..  R  ' 

>  *  ~r  (»E  *  »  -  '  jA 


sin  2y 


(41) 


By  successive  differentiation  these  equations  can  be  combined  to  give  a 
fourth  order  equation  in  7  as  follows 

7  ♦  7  ♦  24kg  9g  sin  27  7  ♦  l8kg  sin  27  -  0  (42) 


While  this  equation  still  requires  that  both  r  and  7  remain  small,  it  does 
not  restrict  the  value  of  7,  Since  it  has  already  been  demonstrated  that 
the  variation  in  7  is  a  low  frequency  solution  with  negligible  damping, 

Eq.  (42)  can  be  simplified  by  dropping  the  first  and  fourth  derivative 
terms.  In  terms  of  the  non-dimensional  time  t  the  equation  becomes 
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from  which  a  first  integral  can  ba  obtained  in  the  fora 


3?  -  -  •iny, 


l  - 


slnS 

HPjTj 


“1 1/2 


(44) 


Equation  (44)  can  then  be  expressed  in  the  fora  of  an  elliptic  Integral  as 
follows 


t  ■  1 


J. 


b 


*L 


I  [>• 


2  2 

sin  yQ  sin% 


P75 


(45) 


where 

6  ■  *1"'1  <!£%>  <“*> 

From  Eq.  (45)  it  is  seen  that  the  resulting  angular  motion  of  the  satellite 
will  be  a  large  amplitude  oscillation  about  the  position  of  the  minor  axis 
of  the  earth' 8  equatorial  section  with  an  amplitude  of  7Q-  The  period  of 
this  osculation  is  given  by 


2 


[l  -  ain2^  sin2*]1/2 


(47) 


where  the  integral  is  K (7Q)>  the  complete  elliptic  integral  of  the  first 
kind. 

The  variations  in  the  orbital  radius  r  during  this  oscillation  can  also 
be  determined  by  substituting  Eq.  (44)  into  Eq.  (41)  and  integrating. 
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Xn  this  Integration  on  the  right  side  of  Eq.  (4l)  is  assumed  to  be 
a 

equal  to  r/\  The  resulting  variation  of  r  is  given  by 

C 


which  can  be  simplified  to 

r  2  f/2 

££  >  siny  1  -  SiBJL  (1*9) 

rc  "  °  L  si?70J 

A  comparison  of  Eqs.  (44)  and  (49)  shove  that 


A£  .  .  I  12 

rc  3  dt 


(50) 


Thus,  when  7  is  increasing  in  magnitude,  Ar  la  negative  and  the  orbital 


radius  is  slightly  less  than  r conversely,  when  7  is  decreasing,  Ar  is 
positive  and  the  orbital  radius  is  slightly  larger  than  r£.  The  maximum 
excursion  in  Ar  occurs  when  7  equals  zero  at  the  minor  axis  and  is  given  by 


♦  4,/kg  sin70 


(51) 
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III.  BESOMS  AMD  DISCU33I0H 

SYNCHRONOUS  SATELLITE 

As  a  result  of  the  foregoing  analysis,  the  effects  of  the  non- spherical 
mass  distribution  of  the  earth  on  a  nominally  synchronous  satellite  can  be 
evaluated.  The  effect  of  the  usual  earth's  oblateness,  or  Jg  term  in  the 
earth's  potential,  is  simply  to  increase  the  radius  of  the  synchronous  orbit 
by  a  small  amount,  as  indicated  in  Eq.  (l6).  *nxe  solution  to  Eq.  (18)  is 
to  a  good  approximation 

rc 

where  r  is  the  solution  of  Eq.  (18)  when  JL  is  zero.  For  ©.,  corresponding 

O  C  Cj 

to  a  period  of  one  sidereal  day  (86l64  sec),  the  value  of  rQ  is  given  by 

r  ■  26194.9  st  ml 
o 

From  Eq.  (52)  it  is  found  that  the  J2  term  Increases  the  orbital  radius  by 
.32  st  ml  to  a  value  of 

rc  -  26195.2  st  mi 

This  effect  would  not  be  difficult  to  eliminate  as  part  of  the  orbital 
injection  corrections. 

On  the  other  hand,  the  effect  of  the  ellipticity  of  the  earth's  equato¬ 
rial  section  is  a  more  serious  problem.  An  examination  of  Eqs.  (37)  and 
(39)  shows  that  for  values  of  A  7  up  to  about  10°  both  equations  can  be  ap¬ 
proximated  by  the  relation 

A7  -  (-  9fc2  ©/sin  27c)t2 


(53) 
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In  a  similar  manner  Eqs.  (36)  and  (38)  can  both  be  approximated  by  the 
relation 

Is  -  ♦  (I2k2  0£  aln  2r0)t  (54) 

c 

From  Eq.  (53)  It  la  possible  to  compute  the  time  required  for  a  satel¬ 
lite  to  drift  from  an  Initial  epoch  angle  7q  through  an  angle  A 7.  Figure  3 
is  a  plot  of  this  time  Interval  for  by  equal  to  10°  as  a  function  of  the 
initial  epoch  angle  yQ.  It  is  seen  that  the  drift  time  has  a  minimum  of 
.1738  yr  or  about  2  months  for  yQ  equal  to  ♦  45°  with  the  drift  being  toward 
the  minor  axis  ( 7  -  0°).  Also,  the  drift  time  becomes  infinite  for  yQ  equal 
to  0°,  +  90°,  l8o°.  However,  ♦  90°  represent  positions  of  unstable  equi- 

O  r  O 

librium  as  demonstrated  by  the  large  amplitude  solution,  while  0  and  loO 
are  positions  of  stable  equilibrium.  Plus,  unless  station-keeping  propulsion 
is  incorporated  into  the  satellite  design  these  two  stable  positions  are 
the  only  locations  at  which  a  truly  synchronous  satellite  can  exist. 

It  is  also  possible  to  determine  by  means  of  Eq.  (54)  the  change  in 
Ar  which  occurs  during  the  10°  drift  which  was  assumed  in  Fig.  3.  In  Fig.  4, 
Ar  is  plotted  as  a  function  of  yQ  and  it  is  seen  that  the  maximum  value  is 

(2) 

about  15  miles.  Thus  the  major  effect  of  the  term  is  to  produce  a 

change  in  the  relative  angular  position,  7,  of  the  satellite  as  seen  from 
the  earth. 

STATION-KEEPING  PROPULSION 

If  it  Is  desirable  to  maintain  a  satellite  in  a  synchronous  position 
other  than  the  two  stable  positions  described  above,  the  amount  of  propulsion 
required  for  station-keeping  becomes  an  important  consideration.  To  determine 
this,  Eqs.  (53)  and  (54)  are  differentiated  to  give 
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Ar  (st  mi) 


Fig.  4  —  Radial  error  as  a  function  of  y^ 
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Lj  m  -  (i8k2  ©*  tin  2y0)t  (55) 

and 

dr  -  -  12k2  ®E  rc  sin  2r0  (56) 

Equation  (56)  gives  tue  radial  velocity  resulting  from  the  perturbation  of 
the  orbit,  while  the  tangential  velocity  component  is  obtained  by  Multiplying 
Eq.  (55)  by  rfi  to  give 

AVe  -  -  l8kg  rc  0*  t  sin  2 jQ  (57) 

From  Eq.  (56)  it  can  be  shown  that  the  constant  radial  component  of  velocity 
has  a  maximum  value  of  .0146  ft/sec  for  yQ  equal  to  45°.  On  the  other  hand, 
from  Eq.  (57)  it  can  be  shown  that  AV^  builds  up  linearly  at  the  rate  of 
•  1395  ft/sec/day.  Thus  after  one  day,  the  tangential  component  of  velocity 
AV0  Is  the  dominant  perturbation  in  the  satellite  mot < on.  To  restore  the 
satellite  to  its  original  position  at  an  angle  jo  after  it  has  drifted  toward 
the  minor  axis  for  time  t  it  is  necessary  to  apply  sufficient  propulsion  to 
change  its  velocity  by  an  amount  2AVQ  in  a  direction  opposite  to  the  drift. 
Since  the  system  is  conservative,  this  reversed  of  the  drift  velocity  will 
cause  the  satellite  to  move  back  to  its  original  value  of  in  an  addition*! 

time  interval,  tQ.  The  satellite  then  reverses  direction  and  begins  to  drift 
again  toward  the  minor  axis  and  the  above  procedure  is  repeated.  Thus,  at 
intervals  of  2tQ,  it  is  necessary  to  make  a  velocity  change  -2AVfi.  In  a 
total  time  of  T  equal  to  2ntQ  the  total  velocity  change  is  given  by 


V,  -  -  2£», 
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“  *  36k2  rc  ©|  ntQ  Bin  2 rQ 

.  ♦  l8kg  rc  ©2  T  8la  2ye  (56) 

where  n  is  the  number  of  corrective  impulses  in  tine  T.  Thus  the  velocity 
change  per  unit  time  is  given  by 

y  m 

r  •  *  l8k2  rc  ®e  •*» 

.  50.9  sin  2 70  (59) 

It  is  seen  that  the  propulsion  requirement  might  be  as  much  as  50.9  ft/ sec 

for  each  year  and  this  figure  is  independent  of  the  time  interval  between 

corrections.  However,  the  time  interval  would  be  determined  by  assigning 

a  permissible  value  of  the  drift  angle  A?  and  using  Eq.  (53)  to  determine 

the  corresponding  value  of  t  . 

o 

It  should  be  emphasised  that  the  above  propulsion  requirement  is  only 
for  station-keeping  and  does  not  include  the  propulsion  required  for  the 
elimination  of  initial  condition  errors  during  orbital  injection. 

A  plot  of  the  yearly  velocity  requirement  superposed  on  a  Mercators 
Projection  of  the  earth  is  shown  in  Fig.  5.  Thus,  for  a  given  position  of 
the  satellite  on  the  equator,  the  yearly  velocity  requirement  con  be  read 
directly.  Based  on  Ref.  1  this  curve  has  been  plotted  with  the  minor  axis 
of  the  earth's  equatorial  section  located  at  123°  9'  West  Longitude,  and 
56°  51 •  East  Longitude. 

The  mass  of  fuel  for  the  above  velocity  requirement  can  be  determined 


from  the  relation 
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pig.  5  —  Velocity  requirement  per  yeor  os  a  function  of  synchronous  position 


2k 


where 

Bp  ■  mass  of  fuel  per  year 
m  -  total  mass  on  orbit 
I  ■  specific  Impulse  of  the  fuel 


(60) 


Thus  for  a  cold  gas  system  (l  *  ?5  sec)  as  much  as  2.1  per  cent  of  the  mass 
on  orbit  would  be  required  per  year  for  station- keeping  propulsion. 

It  should  be  noted  that  in  addition  to  the  fuel  weight,  an  attitude 
sensing  and  control  system  would  be  required  for  this  type  of  operation, 
resulting  in  an  additional  weight  requirement. 


IARQE  AMPLE  OSCILLATIONS 

In  the  absence  of  any  station-keeping  propulsion,  the  satellite  will 
execute  large  angle  oscillations  in  the  equatorial  plane  about  the  position 
of  the  minor  axis  (7  ■  0°  or  l80°).  The  period  of  these  oscillations  is 
given  by  Eq.  (47)  and  has  been  plotted  as  a  function  of  the  amplitude,  , 
in  Pig.  6. 

During  these  oscillations  the  values  of  7  and  dr  as  a  function  of  time 
can  be  determined  from  Eqs.  (45)  and  (49).  Figure  7  is  a  plot  of  &r  as  a 
function  of  7  for  an  oscillation  with  an  initial  amplitude  of  4p°,  with 
time  points  at  intervals  of  one  tenth  of  the  total  period  (Tq  ■  l.J?4l  yr ) . 
It  is  seen  that  starting  at  jQ  equal  to  +  45°  the  satellite  initially  moves 
toward  smaller  values  of  7  and  toward  positive  increments  in  the  orbital 
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radius.  The  increment  dr  reaches  a  maxima  of  ♦  25.9  mi  at  7  equal  to  zero, 
after  vhlch  it  decreases  to  aero  at  -  45°.  As  7  begins  to  increase,  dr 
becomes  negative,  reaching  a  minimum  of  -  25.9  miles  at  7  equal  to  0  and 
then  increasing  to  aero  as  the  oscillation  completes  its  cycle  at  +  45°. 
Figure  8  shows  the  magnitude  of  the  maximum  excursion  of  dr  as  a  function 
of  the  amplitude  70>  as  determined  by  Eq.  (51). 

From  the  above  discussion  it  is  seen  that  a  satellite  in  the  immediate 
vicinity  of  the  minor  axis  (7  >0°  or  l6o°)  will  execute  small  angular  oscil¬ 
lations  about  either  of  these  two  stable  positions.  On  the  other  hand,  a 
satellite  in  the  vicinity  of  the  major  axis  (7*1  90°)  is  in  unstable 
equilibrium  and  the  slightest  disturbance  will  start  a  long  period  oscil¬ 
lation  with  a  double  amplitude  of  l8o°  about  the  position  of  the  minor  axis. 

Reference  2  indicates  that  a  synchronous  satellite  would  be  a  very 

effective  tool  for  the  determination  of  the  ellipticity  of  the  earth's 

equator  al  plane.  It  is  evident  from  the  above  analysis  and  discussion 

that  if  a  satellite  of  this  type  were  established  in  orbit  and  allowed  to 

execute  large  angle  oscillations,  the  midpoint  of  oscillation  would  be  at 

the  position  of  the  minor  axis  of  the  equatorial  section.  A  measurement  of 

the  oscillatory  period  and  the  amplitude,  7  ,  could  be  uaed  in  Eq.  (47)  to 

(?) 

determine  the  value  of  and  thereby  the  parameter  J2'  '  in  the  potential 


function  of  the  earth 
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IV.  CONCLUSIONS 

As  a  result  of  the  analysis  presented  in  this  paper,  the  following 
conclusions  can  be  stated. 

o  The  effect  of  the  earth's  equatorial  bulge  (J2  term  in  the  potential) 
is  to  increase  the  radius  of  a  synchronous  satellite  orbit  by  .32  mile. 

0  The  effect  of  the  ellipticity  of  the  earth's  equatorial  section 
(J2  '  term  in  the  potential)  is  to  produce  large  angle  oscillations  of  the 

satellite  about  the  position  of  the  minor  axis  of  the  equatorial  section 
with  periods  in  excess  of  1.3  years. 

o  Associated  with  these  large  angle  oscillations  are  variations  of  the 
orbited  radius  of  the  same  period  but  with  amplitudes  less  than  37  miles. 

0  Two  stable  positions  exist  for  a  synchronous  satellite.  These  are 
at  the  longitudes  of  the  extremities  of  the  minor  axis  of  the  equatorial 
section,  namely,  123°  9’  W  and  56°  51'  E. 

o  The  longitudes  of  the  two  extremities  of  the  major  axis  of  the 
equatorial  section  are  positions  of  unstable  equilibrium  at  which  any  small 
disturbance  will  set  up  an  angular  oscillation  with  a  double  amplitude  of 
180°  about  the  position  of  the  minor  axis. 

o  To  establish  a  synchronous  satellite  at  any  longitude  other  than  the 
two  stable  positions,  it  is  necessary  to  provide  station-keeping  propulsion 
which  may  amount  to  as  much  as  51  ft/sec/yr.  'Dlls  requirement  is  in  addition 
to  that  necessary  for  the  initial  orbital  injection  corrections. 

o  A  synchronous  satellite  would  be  a  useful  tool  for  determining  a 
more  exact  value  of  the  ellipticity  of  the  earth's  equatorial  section  both 
in  magnitude  and  orientation  of  the  axes. 
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